Spectroscopic Observations and Analysis of the Peculiar SN 1999aa 



G. Garavini 1 , G. Folatelli 1 , A. Goobar 1 , S. Nobili 1 , G. Aldering 2 ' 3 , A. Amadon 4 , R. Amanullah 1 , 
P. Astier 5 , C. Balland 5 ' 6 , G. Blanc 2 - 7 , M. S. Burns 8 , A. Conley 2 - 3 ' 9 , T. Dahlen 10 - 11 , 
S. E. Deustua 212 , R. Ellis 13 , S. Fabbro 14 , X. Fan 15 , B. Frye 2 , E. L. Gates 16 , R. Gibbons 2 , 
G. Goldhaber 2 - 9 , B. Goldman 17 , D. E. Groom 2 , J. Haissinski 18 , D. Hardin 5 , 1. M. Hook 19 , 
D. A. Howell 2 , D. Kasen 2 , S. Kent 20 , A. G. Kim 2 , R. A. Knop 21 , B. C. Lee 2 , C. Lidman 22 , 
J. Mendez 23 - 24 , G. J. Miller 25 ' 26 , M. Moniez 18 , A. Mourao 14 , H. Newberg 27 , P. E. Nugent 2 , 
R. Pain 5 , O. Perdereau 18 , S. Perlmutter 2 , V. Prasad 2 , R. Quimby 2 , J. Raux 5 , N. Regnault 2 , 
J. Rich 4 , G. T. Richards 28 , P. Ruiz-Lapuente 24 , G. Sainton 5 , B. E. Schaefer 29 , K. Schahmaneche 5 , 
E. Smith 21 , A. L. Spadafora 2 , V. Stanishev 1 , N. A. Walton 30 , L. Wang 2 , W. M. Wood-Vasey 2 ' 9 , 
(THE SUPERNOVA COSMOLOGY PROJECT). 



ABSTRACT 

We present an extensive new time-series of spectroscopic data of the peculiar SN 1999aa in NGC 2595. 
Our data set includes 25 optical spectra between -11 and +58 days with respect to B-band maximum 
light, providing an unusually complete time history. The early spectra resemble those of a SN 1991T- 
like object but with a relatively strong Ca H&K absorption feature. The first clear sign of Si n /16355, 
characteristic of Type la supernovae, is found at day -7 and its velocity remains constant up to at least 
the first month after B-band maximum light. The transition to normal-looking spectra is found to occur 
earlier than in SN 199 IT suggesting SN 1999aa as a possible link between SN 1991T-like and Branch- 
normal supernovae. Comparing the observations with synthetic spectra, doubly ionized Fe, Si and Ni 
are identified at early epochs. These are characteristic of SN 1991T-like objects. Furthermore, in the 
day - 1 1 spectrum, evidence is found for an absorption feature which could be identified as high velocity 
C ii /16580 or Ha. At the same epoch C m /14648.8 at photospheric velocity is probably responsible for 
the absorption feature at 4500 A. High velocity Ca is found around maximum light together with Si n and 
Fe ii confined in a narrow velocity window. Implied constraints on supernovae progenitor systems and 
explosion hydrodynamical models are briefly discussed. 

Subject headings: supernovae: general - supernovae: individual (SN 1999aa) 
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1. Introduction 

The observed homogeneity and brightness of Type 
la supernovae (SNe la) make them excellent tools for 
distance estimates over extremely large distances and 
hence for measurements of cosmological parameters 
(see e.g. Perlmutter et al. (1998); Garnavich et al. 
(1998); Schmidt etal. (1998); Riess et al. (1998); Perl- 
mutter et al. (1999); Knop et al. (2003); Tonry et al. 
(2003)). However, cosmological results derived from 
supernovae rely on the evidence that distant explo- 
sions are similar to well-studied nearby ones and that 
they can be calibrated with the same techniques, e.g. 
the luminosity-light curve-timescale relation (Phillips 
1993). The Supernova Cosmology Project (SCP) co- 
ordinated an extensive campaign to study a large num- 
ber of z < 0.1 SNe la in the spring of 1999 in order 
to better understand the intrinsic properties of SNe la 
and thereby improve cosmological measurements us- 
ing them (Aldering 2000; Nugent, Aldering, & The 
Nearby Campaign 2000). The subject of this work, 
SN 1999aa, was one of the SCP targets in that cam- 
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paign. 

At this time, several fundamental questions about 
Type la supernova physics remain. The nature of the 
progenitor system is still poorly constrained as are 
the details of the explosion and thus the origin of 
many of the differences observed among Type la su- 
pernovae. Normal SNe la (sometimes called Branch- 
normal, Branch et al. (1983); Branch, Fisher, & Nu- 
gent (1993)) present early spectra dominated by in- 
termediate mass elements (IMEs), such as Mg, Si, S, 
and Ca, which are replaced by features due to iron- 
peak ions (such as Fe, Co, and Ni) as the spectrum 
evolves with time. However, peculiar events such as 
SN 199 IT, SN 1997br and SN 2000cx (Filippenko 
et al. 1992; Jeffery et al. 1992; Phillips et al. 1992; 
Ruiz-Lapuente et al. 1992) show different character- 
istics. Their early spectra have weak 1ME lines and 
strong doubly ionized iron lines while spectra after 
maximum light look almost completely normal. More- 
over, their light curves are generally characterized by a 
slow post maximum decline rate, and thus, SN 1991T- 
like supernovae are sometimes called peculiar slow de- 
cliners. These characteristics have been regarded as 
possible signs of different classes of progenitors. 

SN 1999aa exhibited spectral characteristics com- 
mon both to Branch-normal and peculiar slow decliner 
SNe la, and thus has been proposed as a key object 
which may help in understanding the physical origin of 
the observed diversity (Branch 2000; Branch, Baron, 
& Jeffery 2001; Li et al. 2001b). 

Hydrodynamical models predict atmospheric com- 
positions of supernovae in good agreement with what 
is found in observed spectra. None of the currently 
available models, though, produces an exhaustive de- 
scription of the whole spectral time evolution and none 
is able to reproduce the complete range of observed 
differences among SNe la. For extensive reviews of 
the theoretical models and observations see Filippenko 
(1997); Hillebrandt & Niemeyer (2000); Leibundgut 
(2000); Livio (2000); Branch et al. (2001). 

The most widely accepted model for Type la super- 
novae involves the thermonuclear disruption of a C+O 
white dwarf star (WD) accreting material from a com- 
panion star (Whelan & Iben 1973; Nomoto 1982; Iben 
& Tutukov 1984; Paczynski 1985). From an observa- 
tional point of view this conclusion is supported by the 
amount of energy released, the absence of hydrogen 
lines (but see Marietta, Burrows, & Fryxell (2000)) 
and the occurrence of SNe la in elliptical galaxies ex- 
clusive of other types. 
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A more controversial issue is whether SNe la are 
the result of explosions at Chandrasekhar mass ( Arnett 
1969; Hansen & Wheeler 1969; Nomoto, Sugimoto, & 
Neo 1976; Khokhlov 1991; Woosley & Weaver 1994) 
or at sub-Chandrasekhar mass (Livne 1990; Livne & 
Glasner 1991; Woosley & Weaver 1994). In the for- 
mer, thermonuclear burning of carbon occurs in prox- 
imity to the center of the star and the burning front 
proceeds outward. In the latter, helium accreting in 
the external layer of the supernova, ignites. A deto- 
nation then propagates outward through the He layer 
and another inward compressing the C+O nucleus that 
ignites off-center. 

While Chandrasekhar mass models have been suc- 
cessful in reproducing many of the observed charac- 
teristics of Branch-normal and SN 1991T-like super- 
novae (Nugent et al. 1997; Fisher et al. 1999), sub- 
Chandrasekhar models are in good agreement with 
fainter explosions such as SN 1991bg-like objects 
(Nugent etal. 1997). 

Pre-maximum spectra can discriminate among 
these two scenarios since sub-Chandrasekhar mod- 
els have external layers dominated by He and Ni and 
do not leave any unburned carbon or produce IMEs at 
expansion velocities above 14000 km s _1 . The identi- 
fication of C, or strong Ca or Si lines in pre-maximum 
spectra would then rule out this possibility. 

Most of the supernovae observed seem to show 
characteristics in agreement with carbon ignition oc- 
curring at the center of the WD. However, this con- 
dition can be the result of both a single degenerate 
C+O WD accreting hydrogen from a companion or a 
merging double degenerate C+O WD, as sometimes 
suggested for SN 1991T-like supernovae (Fisher et al. 
1999). 

The detection of narrow hydrogen lines in super- 
nova spectra (Hamuy et al. 2003) (see also a discus- 
sion in Livio & Riess (2003)) would favor the single 
degenerate scenario but a non-detection cannot rule out 
these models. A substantial amount of hydrogen can 
be removed from the companion star and get mixed 
into the exploding WD (Marietta et al. 2000). In this 
case the detection of H lines is expected even in low 
resolution spectra (Lentz et al. 2002) and would then 
exclude a double degenerate progenitor system. 

While single degenerate models are currently con- 
sidered the most promising, other questions remain 
open. The hydrodynamics of the explosion and the 
details of the flame propagation pose many numeri- 



cal and conceptual problems currently under investiga- 
tion, e.g. computational resolution in 3D simulations, 
and details of flame instabilities, for extensive discus- 
sion see e.g. Hillebrandt & Niemeyer (2000) and Blin- 
nikov & Sorokina (2002). Several models have been 
proposed to describe the explosion mechanism. Next, 
we highlight the possible spectroscopic observables 
that could help in constraining the parameters of the 
models. 

Pure one-dimensional (ID) deflagration models 
(such as W7, Nomoto, Thielemann, & Yokoi (1984)) 
have carbon present down to v > 14900 km s _1 while 
ID delayed-detonation (DD) models (Khokhlov 1991; 
Yamaoka et al. 1992; Woosley & Weaver 1994) burn 
C almost completely up to v ~ 30000 km s" 1 . Thus, 
the identification of carbon lines at low velocity would 
disfavor published DD models. Delayed-detonation 
models can also produce an unusually high Doppler 
blue shift of IMEs by tuning the deflagration to detona- 
tion transition density (Lentz et al. 2001). Thus, IME 
lines confined at higher than normal velocities could 
be more naturally described by DD models. Lines of 
stable Fe and Ni at high velocity in the spectra prior 
to maximum are also consistent with the prediction of 
a deflagration to detonation transition (Iwamoto et al. 
1999). 

The complexity of explosion models has advanced 
to the level where deflagration in 3D can be explored 
(Khokhlov 2000; Gamezo et al. 2003). These results 
motivated a few attempts to investigate the spectral 
outcome of these models by means of parametrized 
synthetic spectra codes (Thomas et al. 2002; Baron, 
Lentz, & Hauschildt 2003; Kasen et al. 2003). In 
3D deflagration, due to the highly convoluted turbulent 
flame propagation, heavy mixing of freshly synthe- 
sized and unburned material can occur (Gamezo et al. 
2003). Evidence for C and O lines at low velocities, as 
well as clumps of burned material (e.g: IME or Fe) in 
the external layers (high velocity), would confirm the 
plausibility of, and the need for, such a degree of com- 
plexity. The same mechanism that mixes carbon and 
oxygen with IME in 3D calculations could probably 
also work for hydrogen. 

Polarization measurements in SNe la spectra (Wang, 
Wheeler, & Hoeflich 1996; Howell et al. 2001; Kasen 
et al. 2003; Wang et al. 2003) have strengthened 
the conviction that some degree of asymmetry can 
be found in supernova atmospheres. In particular, 
the high velocity (HV) component of the Ca n IR 
triplet found in some supernovae (e.g: SN 2000cx, 
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SN 200 lei) required 3D simulations to be fully repro- 
duced (Thomas et al. 2003; Kasen et al. 2003). 

To explore the results of 3D explosion models, full 
three-dimensional radiative transfer calculations of su- 
pernova spectra would be required. However, these are 
not yet available and ID parametrized radiative trans- 
fer calculations are still in their infancy. Thus, direct 
analysis by means of parametrized codes - both in ID 
and in 3D - remains the fastest and most versatile way 
of testing models predictions and guiding new devel- 
opments. 

In the present work we present a new, comprehen- 
sive spectroscopic timeseries for SN 1999aa. We use 
the direct analysis code SYNOW (Jeffery & Branch 
1990; Fisher et al. 1997, 1999) as a tool to describe 
the data and to identify potentially interesting features. 
This, together with the measurements of the velocities 
inferred from the minima of the spectral features, is 
used to investigate the structure of the expanding at- 
mosphere. In particular, we have looked for evidence 
of carbon, oxygen, and hydrogen lines in early spectra 
as well as the velocity ranges of iron, nickel and inter- 
mediate mass element lines, in order to try to answer 
some of the open questions outlined above. 

This paper is organized as follows. The spectro- 
scopic data of SN 1999aa, and a short description of 
the data reduction scheme are presented in section 2. 
The analysis methodology is introduced in section 3. 
In section 4 our SN 1999aa spectra are compared with 
those of spectroscopically peculiar and normal super- 
novae taken from the literature. The synthetic spectra 
for days -11, -1, +5, +14 and +40, produced with 
the highly parametrized SYNOW code, are discussed 
in the same section outlining the spectral peculiarities 
of this object. Velocities inferred from several spec- 
tral features are analyzed in section 5. A discussion 
about possible implications for supernova models is 
presented in section 6 and our conclusions are given 
in section 7. 

2. Data-set and reduction procedure 

SN1999aa in NGC 2595 (R.A. = 8/i27m42".03, 
Decl. = +21°29'14".8, equinox J2000.0) was discov- 
ered independently by Armstrong & Schwartz (1999), 
Qiao et al. (1999) and Nakano, Kushida, & Kushida 
(1999) on February 11, 1999. Filippenko, Li, & 
Leonard (1999) noted that it was a SN 1991T-like Type 
la supernova based on a spectrum taken the follow- 
ing night. Fig. 1 shows the position of SN 1999aa 



in its host galaxy NGC 2595, at redshift z = 0.0146. 
NGC 2595 is an SBc galaxy with a compact blue core 
and blue spiral arms. The recession velocity of the 
host galaxy was determined from narrow H-alpha and 
[N n] emission (Armstrong & Schwartz 1999). Ac- 
cording to Schlegel, Finkbeiner, & Davis (1998) the 
Galactic reddening in the direction of SN 1999aa is 
E(B - V) = 0.04 mag. Based on the SN colors, Krisci- 
unas et al. (2000) concluded that host galaxy redden- 
ing is negligible. They also estimated the light curve 
decline rate Am^ =0.746+0.024 from a fourth-order 
polynomial fit to the light curve. Jha (2002) reported 
Am 15 =0.85±0.08 using a different photometry dataset. 

The SCP follow-up campaign of this supernova in- 
volved 5 different instruments and resulted in 25 op- 
tical spectra ranging between 1 1 days before and 58 
days after maximum light (all epochs in this work are 
given with respect to the B-band maximum). In most 
cases the spectra were acquired using different instru- 
mental settings for the blue and the red parts of the 
spectrum in order to avoid possible second order con- 
tamination. Whenever both blue and red spectra were 
taken at the same epoch, we present the combined 
spectrum. The observations were performed align- 
ing the slit along parallactic angle in order to min- 
imize light loss from differential atmospheric refrac- 
tion. Specifications of the data-set are provided in Ta- 
ble 1 and the spectral time sequence is shown in Fig. 2. 

All the raw data were analyzed with a common re- 
duction scheme using standard IRAF routines. The 
two-dimensional spectra were bias-subtracted and flat- 
fielded using calibration images taken with the same 
instrument settings as the SN spectra and during the 
same night of observation. In most cases, the obser- 
vations were split into multiple exposures in order to 
allow elimination of cosmic rays from the final spec- 
trum. Background subtraction was performed on the 
resulting spectra using a fitted model of the underlying 
sky and host galaxy signal. The supernova spectrum 
was extracted using the variance weighted optimal 
aperture extraction method (Home 1986). Several arc 
lamp exposures, generally taken off supernova posi- 
tion, were used for wavelength calibration each night. 
The accuracy of the calibration was checked against 
sky lines and no significant deviations were found. 
Atmospheric extinction corrections were applied us- 
ing tabulated average extinction curves provided by 
the observatories. Spectrophotometric standard stars 
(Oke 1990; Massey & Gronwall 1990; Turnsheket al. 
1990; Hamuy et al. 1992, 1994; Bohlin et al. 1995; 
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Stone 1996; Bohlin et al. 2001) were observed dur- 
ing each night and their latest calibrations 1 were used 
to flux calibrate the SN spectra. A comparison be- 
tween observations of different standard stars during 
the same night was used to check for possible system- 
atic errors. A correction for Galactic extinction was 
performed using the standard procedure in Cardelli, 
Clayton, & Mathis (1989) assuming R v = 3.1. Finally, 
residual host-galaxy contamination was checked sta- 
tistically. The spectrum of the host galaxy at the posi- 
tion of the supernova, multiplied by an arbitrary scale 
factor, was subtracted to the observed data and com- 
pared to a supernova spectral template. The value of 
the scale factor was determined by minimizing the x 2 
between the model (the supernova spectral template) 
and the data (supernova data plus rescaled host galaxy 
data). The host galaxy contamination was found to be 
negligible at all studied epochs. Thus no (additional) 
galaxy spectrum was subtracted. Telluric absorptions 
and residual fringing patterns were left uncorrected be- 
cause they do not affect our present analysis. For a 
complete description of the data reduction methodol- 
ogy see Folatelli (2003). 

3. Spectral Analysis Methodology 

As we will use SYNOW to produce synthetic spec- 
tra and to investigate line identifications and velocity 
ranges of ions in SN1999aa, here we briefly review 
its underlying precepts and parameters. SYNOW gen- 
erates spectra using a simple conceptual model of an 
homologously expanding supernova envelope. This 
model consists of a continuum-emitting, sharply de- 
fined photosphere surrounded by an extended line- 
forming, pure scattering atmosphere. Line transfer is 
treated using the Sobolev method (Sobolev 1960; Cas- 
tor 1970; Jeffery & Branch 1990). Thus, line opac- 
ity is parametrized in terms of Sobolev optical depth. 
The choice of ions used in the calculation is guided by 
experience and the SN ion signatures atlas of Hatano 
et al. (1999b). For each ion introduced, Sobolev opti- 
cal depth as a function of radius for a "reference line" 
(usually a strong optical line) is specified. Optical 
depths in other lines of the ion are set assuming Boltz- 
mann excitation of the levels at temperature T exc . 

The parameters \ p hot an d Tbb set the velocity and 
blackbody continuum temperature of the photosphere, 
respectively. For each ion, the optical depth r and 

'We have corrected the (Hamuy et al. 1992, 1994) spectra for telluric 
atmospheric features. 



the specified minimum ejection velocity v ml „ is given. 
The optical depth scales exponentially with velocity 
according to e-folding velocity, v e , up to a maximum 
velocity given by w max . If w min >w pho , for an ion, we 
refer to the ion as "detached." 

The black body assumption is a basic simplification 
of the processes that contribute to form the continuum 
emission, therefore, r bb cannot be regarded as physi- 
cal information. Thus, SYNOW produces only a rough 
indication of the continuum level, which can disagree 
with the observed spectrum. We handled this by bring- 
ing the blue portions of the observed and modeled 
spectra into agreement. The continuum mismatches 
which then occur at longer wavelengths have not ad- 
versely affected our study of the few spectral features 
present in the red. 

Direct analysis results, such as those from SYNOW, 
yield meaningful constraints for hydrodynamical ex- 
plosion modelers. Furthermore, direct analysis often 
reveals the presence of lines that otherwise go unde- 
tected without some treatment of the line blending that 
characterizes supernova spectra. 

4. Spectral Data and Modeling 

Prior to B-band maximum light, normal SNe la 
have spectra dominated by intermediate mass elements 
such as Si n, S n, Mg n, Ca n and O i (Branch et al. 
1983, 1993). With time, the absorption features from 
these elements become weaker and increasingly con- 
taminated by iron-peak element lines. The substitu- 
tion is usually complete around 30 days after maxi- 
mum, when the photosphere of the supernova starts 
receding into the iron-peak core (i.e. the innermost 
part of the ejecta where mainly iron-peak elements are 
present). Objects such as SN 1991T, SN 1997br and 
SN 2000cx represent a deviation from the impressive 
homogeneity of the spectral and photometric charac- 
teristics of SN explosions (Fisher et al. 1999; Li et al. 
1999, 2001a). This has raised the question whether 
they could be explained as a different physical phe- 
nomenon or as extreme cases of the same process. 
Ever since the first spectra of SN 1999aa started cir- 
culating in the supernova community, several authors 
suggested that this object could be helpful in address- 
ing this issue (Branch 2001; Li et al. 2001b). This 
work should be viewed in that context. 

In order to understand and interpret the differences 
between SN 1999aa and other normal and peculiar su- 
pernovae we now analyze the spectral time evolution 
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through spectral comparison and spectral modeling. 
We have selected five epochs, approximately one week 
apart from each other, to describe the spectral time 
evolution in all interesting phases. These are: -11,-1, 
+5, +14, and +40 days with respect to B-band maxi- 
mum. 

4.1. Day 11 before maximum 

Our spectrum of SN 1999aa at 1 1 days prior to B- 
band maximum light is shown in Fig. 3. The spec- 
tra of SN 199 IT, SN 1990N and SN 1994D are also 
shown for comparison. The identification of the lines 
labeled here and in the following graphs follow those 
in Li et al. (1999, 2001a); Fisher et al. (1999); Patat 
et al. (1996); Mazzali, Danziger, & Turatto (1995); 
Kirshner et al. (1993); Jeffery et al. (1992). The spec- 
tra of SN 1 999aa and SN 1 99 1 T are very different from 
those of normal supernovae SN 1990N and SN 1994D 
around 1 1 days prior to the B-band maximum light. In- 
stead of the typical Si n, S n and Mg n lines, early spec- 
tra of these peculiar SNe la are dominated by two deep 
absorptions due to Fe m AA4404 and 5 129. Si m i4560 
is responsible for the smooth absorption on the red 
side of the Fe in /14404 line. Possibly weaker Ni m 
lines contribute near 4700 A and 5300 A, (Jeffery et al. 
1992; Ruiz-Lapuente et al. 1992). What mainly distin- 
guishes SN 1999aa from SN 199 IT is the presence in 
the former of the trough around 3800 A, most probably 
due to Ca n H&K. This absorption is weak compared 
to normal SNe la. The weak and broad line around 
6150 A could be due to Si n /16355 with some contam- 
ination from C ii /16580. 

Synthetic spectrum 

In Fig. 4 we show the results for the best match 
synthetic spectrum computed using the parameters 
in Table 2. The continuum black body temperature 
is set to 13700 K and the photosphere is placed at 
11000 km s _1 . The dominant and characteristic fea- 
tures at this epoch are due to Fe in lines. While some 
authors identify the small absorption line at 4000 A 
only as Si n /14130, in our synthetic spectrum we find 
that Si ii alone cannot reproduce the feature completely 
due to its small optical depth. If we add a contribution 
from Co ii (as in Hatano et al. (1999a) for example) 
we are able to reproduce the observed spectrum. Note 
that to improve the match Co n needs to be detached 
above the photosphere. Si m appears confined below 
16 x 10 3 km s _1 (similarly to Si n) and is responsible 



for the features at 4400 A, and 5500 A, as originally 
identified by Jeffery et al. (1992) for SN 1991T and 
SN 1990N. 

The only way we could account for the complete 
structure of the Ca ii H&K features and for its velocity 
distribution (blue edge at 40000 km s~') was to intro- 
duce a HV component (with higher v e ) in addition to 
the one at v max = 19500 km s _1 . This allows us to 'fill- 
up' the whole spectral profile at this epoch and at day 
-1 it will produce the observed blue absorption feature 
of the Ca n IR triplet, as noticed already in Hatano et al. 
(1999a) for SN 1994D and analyzed in Kasen et al. 
(2003) for SN 2001el and in Thomas et al. (2003) for 
SN 2000cx. The region around 4600 A and the broad 
absorption at 5300 A can be matched with detached 
Ni in. 

We have focused in on the wavelength regions 
around 6150 A and 4500 A for particular study, since 
these have been discussed for other SNe la in the liter- 
ature. Si ii 6150 has a very broad profile. As we will 
show in section 5 the velocity fitted for Si n 6355 at 
this epoch is the lowest among all known supernovae 
(v~8000 km s~'). This suggests that probably another 
ion is partially or totally responsible for this absorp- 
tion feature. In Fig. 3, a similar absorption feature 
is visible for the spectrum at -10 days of SN 199 IT 
Fisher et al. (1999) discussed the possible presence of 
HV C in early spectra of SN 199 IT based on the line 
at 6150 A usually attributed to Si n 6355. The same 
possibility has been also proposed for SN 1990N in 
Fisher et al. (1997) and Mazzali (2001), for SN 1994D 
in Hatano et al. (1999a) and recently for SN 1998aq in 
Branch et al. (2003). 

Guided by the low velocity value of Si n 6355 and 
by the similarity with SN 199 IT we also introduced 
a HV C ii (v m ,„=19 x 10 3 km s~') in order to try to 
improve the description of the line profile at 6150 A. 
Fig. 5 shows different S YNOW spectra compared with 
the data in the 6150 A region. The small notch around 
6200 A is an artifact of the spectrum extraction and 
does not have to be considered a supernova feature. 
The error spectrum on the wavelength region 6195- 
6210 A is 30% larger than that of neighboring regions 
thus, the small notch does not have large statistical 
weight for any quantitative measurements. The contin- 
uum level of the models has been shifted to match the 
data in this region. The model at the top has no C and 
includes only Si n. The absence of a C n component, 
i.e. considering only Si n, places the minimum of the 
feature to the blue side of the absorption at 6150 A and 
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does not reproduce completely the line profile. The pa- 
rameter V m in is set equal to the photospheric velocity 
and thus has been determined looking at all the lines 
together. The dashed line indicates the synthetic spec- 
trum with the maximum velocity of Si n set to a value 
(Vmax = 30xl0 3 km s -1 ) higherthan in the later epochs, 
but which reproduces more accurately the line profile. 
The model at the bottom of Fig. 5 includes both car- 
bon and silicon. Because of the introduction of C n, we 
would expect a small absorption (C n /17234) around 
6800 A (close to a telluric line) which is not convinc- 
ingly visible. It is not fully evident which of the two 
models proposed matches better the data. The result 
of a% 2 analysis shows that they are equally plausible. 
However, the introduction of C n would solve the puz- 
zle of the low velocity value measured for Si n -16355. 
Note that later on (at day -7, Fig. 2) the shape of the 
feature at 6150 A changes with the appearance of the 
first strong Si n 16355 line. 

In our spectrum of SN 1999aa we find that the ab- 
sorption on the red of Si m 14568 can be well matched 
by C in /14648.8 confined to low velocities. Its con- 
tribution is shown in Fig. 6 where synthetic spectra 
with and without a C in component are compared with 
the observed spectrum. This absorption is unusual for 
SNe la and its identification is still unclear. Three 
dimensional models indicate that unburned doubly- 
ionized C can, in principle, be visible in the hot interior 
of the supernova atmosphere (Gamezo et al. 2003). 
Later in the spectral evolution, the identification of C 
could be difficult due to increased blending of spectral 
lines. Further, we note that this absorption feature is 
no longer visible in our spectrum at day -7 (see Fig. 
2). 

As an alternative to C in and C n, H could be used 
to match the absorption at 4500 A and 6150 A (Lentz 
et al. 2002). We discuss these alternative explanations 
in section 6.3. 

4.2. Day 1 before maximum 

The spectrum of SN 1999aa around maximum light 
is shown in Fig. 7 together with those of SN 199 IT, 
SN 1981B and SN 1994D. SN 1999aa resembles that 
of normal SNe la, while SN 199 IT is still dominated 
by Fe m lines with very weak Si n and Ca H&K lines 
at this epoch. The Ca n H&K region of SN 1999aa 
shows a characteristic split where the red component 
is gaining strength. S n appears at this epoch show- 
ing the typical 'W shaped feature at 115454,5606 A 



that SN 1991T-like SNe la do not usually show. 
Si ii 16355 is visible but fainter and redder than in 
normal SNe la. If the evolution from C n to Si n pro- 
posed by Fisher et al. (1999) is confirmed, SN 199 IT 
at this epoch shows a weaker Si n 16355 absorption 
than SN 1999aa. This suggests that the amount of Si n 
present in SN 1999aa and/or the atmospheric tempera- 
ture could be between that of normals and SN 1991T- 
like SNe la. The Ca n IR triplet around 8000 A has 
two weak features not present in SN 198 IB but visible 
in SN 1994D. All the other major features are the same 
as those seen in normal SNe la. 

Synthetic spectrum 

The synthetic spectrum for day -1 is shown Fig. 
8. As in the previous epoch, Si n and Si in have been 
introduced with a low maximum velocity in order to 
match the observed wavelengths of the minima; see 
Table 3 for details. The S n, Mg n, Co n and Fe n lines 
are well reproduced in detached layers. 

The Ca H&K feature is formed by a blend of the 
two dominant components of Ca n and weaker contri- 
butions from Si n, Co n and Ni n. The double compo- 
nent of Ca ii is required at this epoch both in the H&K 
and IR triplet regions (respectively left and right panel 
in Fig 9). The split seen at the location of Ca H&K 
has been reproduced with Si n 13858 or Si m -13801 
in NLTE simulations (Nugent et al. 1997; Lentz et al. 
2000) but with SYNOW, the only way to account for 
the blue minimum of Ca H&K is to add a high ve- 
locity component (Ca n HV) (Hatano et al. 1999a). 
Ca ii HV with v m ; n = 21000 km s _1 also matches the 
two small absorptions in the blue part of the Ca ii IR 
triplet. In the right panel of Fig. 9 the black body con- 
tinuum has been shifted to the level of the observed 
spectrum around 8000 A to show the match in veloc- 
ity. The correct relative intensity of the two HV Ca ii 
IR triplet features and the correct line profile are better 
reproduced by modeling the 3D geometry of the HV 
clump (Kasen et al. 2003; Thomas et al. 2003). In this 
work we are only interested in identifying the velocity 
range at which this component is present. 

The absorption features we tentatively identified as 
C ii and C m in the early spectrum are no longer de- 
tectable. The broader and stronger Si n 16355 and 
Si in -1505 1 now dominate the wavelength ranges at 
which these lines were present. 

As seen in Fig. 8, the synthetic spectrum fails to de- 
scribe the data in several wavelength regions: around 
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4600 A, 5500 A and beyond 6200 A. In similar analy- 
ses with SYNOW, attempts were made to improve the 
match around 4600 A by adding a high-velocity com- 
ponent of Fe n. For SN 1999aa, this approach intro- 
duced mismatches in other wavelength regions. The 
disagreement, as well as that around 5500 A, could 
also be (partly) explained by the inaccuracy of the 
blackbody continuum, which is clearly showing a mis- 
match for the flux redward of 6200 A. 

As for the spectrum at day -11, the black body 
continuum fails to produce the correct flux redward of 
6200 A but this does not affect the reliability of the line 
profile reproduction. 

Despite attempting several different parameter 
combinations, the synthetic spectrum shows signifi- 
cant disagreement with the data at 4600 A on the red 
side of S ii. 

4.3. Day 5 after maximum 

The spectrum of SN 1999aa around one week af- 
ter maximum, is shown in Fig. 10 and compared with 
those of SN 1991T, SN 1994D and SN 1990N. The 
spectrum of SN 1999aa is quite normal, now also 
showing the blue component of the Si n duo in the 
5700-6200 A region. SN 199 IT still exhibits weaker 
Si ii and S n lines at this epoch. In normal SNe, ab- 
sorption features due to elements such as Co n, Fe n, 
Ca ii IR triplet and O i A1113 get stronger at this time. 
From Fig. 10 and from the spectra of other SN 1991T- 
like supernovae shown in Li et al. (1999, 2001a), it 
is clear that the general trend for these objects at this 
epoch is to still have very weak intermediate mass el- 
ement components and a persistent dominant presence 
of doubly ionized Fe. Furthermore, SN 1997br and 
SN 2000cx (Li et al. 1999, 2001a) already have a large 
contribution from low ionization iron-peak elements 
(more prominent than in Branch-normal SNe), with 
SN 1997br more advanced in this transition. This sug- 
gests a very interesting sequence among these objects 
in which the duration in the phase of the intermediate 
mass elements is shortest for SN 1997br, and longest 
(and thus more similar to normal SNe) for SN 1999aa, 
with SN 199 IT placed somewhere in the middle. For 
Branch-normal supernovae, this phase extends from 
well before maximum to a couple of weeks after. In 
the case of these peculiar objects it seems instead to 
vary in duration and phase during the time evolution. 
This is probably due to the ratio of Fe-peak to IME 
abundances, their radial distribution and temperature 



profile. Further studies are needed to quantify these 
differences. 

Synthetic spectrum 

In the synthetic spectrum for day 5 after maximum, 
shown in Fig. 11, doubly ionized Si and Fe have lower 
optical depth than in the previous epoch analyzed; see 
Table 4. Fe m and Si m are now blended, forming a 
deep trough at 4400 A together with Co n and Mg n. 
We also note a mismatch between the data and the syn- 
thetic spectrum around 4600 A. 

In order to reproduce the three distinct minima be- 
low 5000 A, the Fe n layer has to be detached, with a 
similar velocity range as in day -1, otherwise the lines 
tend to blend together forming a single broad feature. 

At this epoch, Si n and Si m also become detached, 
as did S n and Mg n in the previous epoch. Now 
they all appear confined above 10 x 10 3 km s -1 . This 
is a probable sign that complete silicon burning has 
stopped at this velocity. 

The feature at 5700 A is a blend of Na i with a 
weak Si n component. The Si n component is weak be- 
cause of the low excitation temperature used; a higher 
temperature would create a redder feature than the one 
seen in our data. 

Unfortunately at this epoch our spectrum does not 
cover the region below 3900 A and the Ca IR triplet is 
too weak to fully determine whether there is still need 
for a HV component. 

4.4. Day 14 after maximum 

Our spectrum of SN 1999aa at day 14 is shown in 
Fig. 12 together with those of SN 1991T, SN 1994D 
and SN 1990N. SN 1999aa is now completely in- 
distinguishable from a normal supernova, such as 
SN 1990N, with the possible exception of a slightly 
weaker Si n /16355. This is to be expected, since small 
differences are visible even among so called normals. 
In comparison SN 199 IT is still completing its inter- 
mediate mass elements phase; S n has not yet blended 
with Na i as in normal objects. The overall contribu- 
tion of S ii, Si ii and Ca n is consistently weaker than 
in typical type la. 

Synthetic spectrum 

In the synthetic spectrum at day +14, shown in Fig. 
13, most of the lines can be identified as Fe ii (Table 
5). As in the previous synthetic spectrum, in order to 
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reproduce the three different minima in the Fe n blend 
region (4800 A) we use a minimum velocity higher 
than that of the photosphere (Table 5). 

The deep line at 5700 A is dominated by Na i, and 
its shape can be reproduced if Na i is introduced with 
higher v e than that used for the other ions. 

During the entire spectral evolution, the flux level 
and shape in the region between 6500 A and 8000 A 
has become increasingly more difficult to match. At 
this epoch, the discrepancy between data and synthetic 
spectrum is fairly evident. This suggests an increasing 
inaccuracy of the assumption of an underlying black 
body continuum. 

4.5. Early Nebular Phase Spectra 

The day +25 days of SN 1999aa is shown in 
Fig. 14. For comparison SN 1991T, SN 1994D 
and SN 198 IB are also shown. The main differ- 
ences in the observed spectra lie in the region around 
6000 A, and are probably due to different Si abun- 
dances among the different SNe la. SN 199 IT shows 
the Fe n /L16238, 6246,6456,65 16 lines and just a very 
weak Si n line in the central part of that wavelength 
region. In SN 1999aa, Si n is more evident, and it 
becomes even stronger in SN 1994D and SN 198 IB. 
The trough near 5700 A due to Si n i5972 and Na i D 
is strong in all the SNe. At the red end of the spectra, 
the four SNe show the typical Ca n IR triplet as a very 
deep and broad absorption. SN 1999aa shows a sec- 
ond minimum in the red part of this feature not firmly 
identified, but possibly consistent with an absorption 
feature due to Co n. 

The spectrum of SN 1999aa at day +40 is shown 
in Fig. 15 and compared with those of SN 199 IT, 
SN 1994D and SN 198 IB. At this stage the spectra are 
all dominated by Fe n and Co n lines formed mainly in 
the deep layers of the atmosphere. Only small differ- 
ences are now visible in the depths of some lines. 

Synthetic spectrum 

At day +40 the supernova has started entering the 
so called nebular phase and SYNOW's assumption of 
a sharp photosphere becomes less physically realis- 
tic. Even though SYNOW was designed only to re- 
produce spectra in the photospheric phase, several au- 
thors have drawn conclusions based on its matches to 
nebular spectra (see e.g. Fisher et al. (1999)). In the 
interest of comparison with their work, we also present 
the synthetic spectrum for this epoch. 



With the photosphere now into the iron-peak core, 
at this epoch the supernova spectrum is mainly formed 
by Fe n and Co n; see Fig 16. The only signs of IMEs 
are Ca n H&K and the Ca IR triplet. As proposed by 
Fisher et al. (1999) for SN 199 IT, the strong and wide 
feature near 7000 A can be reproduced by forbidden 
lines of [O n] OM7320,7330) with a high v e . Note that 
in SYNOW this line is also treated with a resonance 
scattering source function. For the details of the com- 
position, see Table 6. Special care was taken in the 
way the components of Fe n and Co n have been in- 
troduced. At this epoch, the spectrum is formed in the 
Fe-peak core but still a part of the outer atmosphere 
makes a contribution. The optical depth of these ele- 
ments in the two regions can be different. 

Mazzali et al. (1998) demonstrated that the expan- 
sion velocity of the Fe-peak core and the luminosity of 
SNe la are correlated. Slow decliner supernovae, such 
as SN 199 IT, are expected to have a higher Fe-peak 
core limit velocity and thus a larger region dominated 
by nuclear statistical equilibrium. In this respect, it 
is interesting to attempt the identification of the iron- 
peak core limit by means of SYNOW synthetic spec- 
tra. Fisher et al. (1999) and Hatano et al. (2002) intro- 
duced an optical depth discontinuity for Co n and Fe n 
at 10000 km s" 1 and 12500 km s" 1 for their synthetic 
spectra ofSN 199 IT and SN 1997br respectively. The 
velocity at which this change occurs can be thought of 
as the iron peak core limit. We tried to reproduce this 
for SN 1999aa. The best match was achieved with a 
discontinuity at 10000 km s -1 . 

5. Expansion Velocities 

The expansion velocities as computed from fits to 
the minima of the spectral lines can provide help in 
investigating the physics of the supernova explosion. 
Weighted fits of the observed minima were performed 
using a non-linear Marquardt-Levenberg minimization 
procedure (Marquardt 1963) applied to a Gaussian 
profile model. The minimum of the line was consid- 
ered to be the center of the Gaussian and the fit un- 
certainty its statistical error. This statistical error is 
usually of the order of a few km s _1 , so is not shown 
on the graphs. In general, the fit was performed on 
the entire absorption of the P-Cygni profile when this 
was well reproduced by a Gaussian model, e.g. when 
there was a low level of contamination by other lines. 
In the contaminated cases the fit was done on a sam- 
ple closer to the bottom of the line. In Figures 17-20 
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we show the time evolution of the velocities for Ca n 
H&K, Si ii ,16355, and Fe m ,U4404,5129 lines for 
SN 1999aa along with several other supernovae from 
the literature. 

Ca ii H&K velocities for SN 1999aa, shown in Fig. 
17, are in the range of those of normal SNe. This fea- 
ture shows an evident split in the minimum at day -3 
and -1. In these cases the measurements have been 
performed selecting a wavelength range that covers 
the whole line profile from approximately 3600 A to 
3850 A. This was chosen in order to obtain a value 
of the velocity that could be compared with the other 
supernovae. Measuring the two minima and deriving 
from those the velocities of the two components of 
Ca ii used in the synthetic spectrum would not be ap- 
propriate since the position of the troughs are the result 
of the blending of Ca n, Si n, Co n and Ni n. 

Si ii velocities are shown in Fig. 18. If the ab- 
sorption near 6150 A of SN 1999aa was due only to 
Si n /16355, the velocities would have been monoton- 
ically decreasing with time, as is seen for the nor- 
mal SN 1994D, SN 1992A and the sub-luminous 
SN 1999by and SN 1991bg. For SN 1999aa the first 
point, 1 1 days before maximum light, has the lowest 
velocity. This is consistent with another ion (proba- 
bly C ii or H according to the SYNOW fit, see sec- 
tions 4.1 and 6.3) also being responsible for the ab- 
sorption feature at this epoch. The wavelength of the 
line minimum remains practically constant during the 
first 20 days after maximum. This is usually inter- 
preted as if the element layer in the supernova atmo- 
sphere is confined to a region above the photosphere, 
in this case around 10100 km s _1 as confirmed in the 
SYNOW synthetic spectra. The same trend is shown 
by SN 2000cx and SN 199 IT, with velocities around 
12000 km s _1 and 9400 km s _1 , respectively, as indi- 
cated by the dashed horizontal lines in Fig. 18. 

The velocities measured for the two Fe m fea- 
tures are shown in Fig. 19 and 20. These lines 
are characteristic of the pre-maximum spectra of the 
SN 1991T-like SNe la and disappear within the first 
week after maximum light. In both cases SN 1999aa 
shows lower velocities compared to SN 199 IT, SN 
1997br and SN 2000cx, although all have a sim- 
ilar slope. Fe in /14404 velocities (Fig. 19), are 
smaller by 2500 km s" 1 with respect to SN 199 IT, and 
3500 km s~ l with respect to SN 2000cx. Fe in /15129 
velocities (Fig. 20) are still smaller, but much closer 
to the velocities of the other peculiar supernovae. 
This could be due to a higher contamination from 



Si ni5051 orFeni5018. 

Generally, the velocity ranges of the lines of 
SN 1999aa are consistent with those of normal su- 
pernovae. The major peculiarity is the restricted atmo- 
sphere region where Si n is present which makes this 
ion appear clearly only one week before maximum 
light. 

6. Discussion 

The analysis carried out in section 4 indicates 
that SN 1999aa was a peculiar object prior to maxi- 
mum light, developing toward normal-looking spec- 
tra around maximum. The SYNOW synthetic spectra 
have shown the velocity ranges in which each adopted 
ion appears, and offer a model of the structure of the 
expanding atmosphere. In this section we summarize 
our findings for different ions and their velocity distri- 
butions and discuss the consequences for constraining 
the progenitor system and the explosion models. 

6.1. Intermediate mass elements 

The velocity range in which IMEs are found con- 
strains the atmospheric region in which incomplete Si 
burning takes place, and thus the hydrodynamics of the 
supernova explosion. We identify the contributions of 
ions such as Si n, S n, Ca n, Na i and Mg n which 
are common to both normal and SN 1991T-like super- 
novae. However, the optical depths and velocity ranges 
found for these elements for SN 1999aa differ from 
those typically used for Branch-normal SNe. The pres- 
ence of unburned C and IMEs (HV Ca n and possible 
signs of C ii) above 14000 km s _1 , as found in our syn- 
thetic spectra of SN 1999aa, suggests that this super- 
nova was not the result of a sub-Chandrasekhar mass 
explosion, as the composition resulting from the mod- 
eling of such explosions do not include any IME or C 
in the external layer of the envelope (see e.g., Livne & 
Arnett (1995) and Woosley & Weaver (1994)). 

The time evolution of the Doppler blue-shift of 
Si ii, Fig. 18, around and after B-band maximum 
shows that this ion is confined within a narrow velocity 
range (v min -10100 km s" 1 ), not only for SN 1999aa 
but also in SN 2000cx (v min -12000 km s _1 ) and 
SN 1991T (v min -9400 km s" 1 ). Our SYNOW mod- 
els of SN 1999aa provide independent evidence for a 
confined Si n layer. Velocities between 10000 km s _1 
and 15000 km s _1 were needed to match the observa- 
tions when the photospheric velocity dropped below 
9500 km s _1 , (see Tables 2 to 6). The same holds for 
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SN 1997br (Hatano et al. 2002) and for SN 199 IT 
(Fisher et al. 1999). For a normal Type la supernova, 
such as SN 1994D, the Si-rich layer extends from the 
photosphere out to v > 25000 km s _1 , (Hatano et al. 
1999a). For SN 1991T-like or SN 1999aa-like ob- 
jects, the weakness of the IME's absorption and the 
earlier than normal domination of the Fe-group ele- 
ments suggest that the region of incomplete Si-burning 
could be confined to a small velocity window. Thus, a 
higher deflagration to detonation transition density in 
DD models could explain some of the spectral pecu- 
liarities. 

In delayed detonation hydrodynamical models this 
would imply an over-luminous object. We do not have 
a sufficiently accurate absolute magnitude measure- 
ment of SN 1999aa, but if the light curve-width lumi- 
nosity relation holds for this object, SN 1999aa should 
be more luminous than normal SNe. Based on the 
measurement of the distance to SN 1991T, Saha et al. 
(2001) and Gibson & Stetson (2001) showed that that 
supernova is not necessarily overluminous. Thus, the 
possible correlations between the absolute magnitude 
of SNe la and the velocity range in which IMEs are 
present requires further studies. 

6.2. Nickel and other doubly ionized elements 

The ionization level of each component of the at- 
mosphere in type la SNe is an indicator of the energy 
balance during the supernovae explosion. Both ther- 
mal and non-thermal ionization are known to be im- 
portant in type la SNe (see e.g. Lucy (1991); Nugent 
et al. (1995); Baron et al. (1996)). 

We have found evidence for Si m, Fe m and Ni m 
in the early spectra of SN 1999aa. As discussed in 
Nugent et al. (1995), doubly ionized Fe and Si in 
the early spectra suggest higher temperature than in 
Branch-normal SNe. The absorption features from 
these ions disappear after day 5, as the supernova at- 
mosphere cools. 

The innermost layers of deflagration models, such 
as W7, are considered to have gone through com- 
plete silicon burning, leaving only a nuclear statistical 
equilibrium composition such as that shown in Hatano 
et al. (1999b). The decay of the resulting 56 Ni to 56 Co 
and then to stable 56 Fe would imply the presence of 
Co m lines not found in the present analysis. As in 
Hatano et al. (2002), we suggest that the Ni m lines 
in our pre-maximum spectra are produced by 54 Fe and 
58 Ni synthesized during incomplete and complete sil- 



icon burning. If this is the case, it is not necessary 
to produce a model that extends the complete silicon 
burning into regions of intermediate mass elements or 
above, as originally proposed for SN 1991T-like ob- 
jects. Among the current DD models, none is able to 
extend the presence of these ions sufficiently far out 
(v max = 30000 km s _1 , Table 2) to match our observa- 
tions. 

Furthermore, the measurements of the Doppler shift 
of Fe in, Fig. 19 and 20, show a lower value than other 
SN 1991T-like objects. This suggests that Fe in ex- 
tends further out in SN 1999aa than in spectroscopi- 
cally normal SNe la, but has its maximal absorption at 
lower radii than other SN 1991T-like objects. 

Hoeflich, Wheeler, & Thielemann (1998) showed 
that a higher initial metallicity would increase the 
abundance of 54 Fe and 58 Ni, and this could possibly 
explain our observations. UV photometry and spec- 
troscopy would yield important information to settle 
this issue and aid in the development of new models 
that could better match the observations (Lentz et al. 
2000). 

6.3. Carbon n or Hydrogen? 

The identification of H or C would impose con- 
straints on either the nature of the progenitor system 
or the hydrodynamics of the explosion respectively. In 
section 4. 1 we suggested the presence of a HV C (C n) 
component for the spectrum of SN 1999aa at day -11. 
In what follows, we discuss possible alternative expla- 
nations, assuming an additional ion is necessary to im- 
prove the matching. 

The evidence for an external layer of carbon in 
Type la supernovae has been discussed in the case 
of SN 199 IT by Fisher et al. (1997) and Fisher 
et al. (1999), for SN 1990N by Mazzali (2001), for 
SN 1994D by Hatano et al. (1999a) and recently for 
SN 1998aq by Branch et al. (2003), based mainly on 
the possible presence of a C n absorption on the red 
side of Si n i6355. For SN 1999aa, this ion helps in 
reproducing the line shape and would explain the ve- 
locity time evolution of the feature at 6150 A, as shown 
in Figs. 4 and 5, and Fig. 18 (section 5). However, due 
to the small optical depth for C n (t = 0.01, table 2) the 
C n 7234 line, which could confirm the presence of C, 
is too weak in the synthetic spectrum to be positively 
identified in the spectrum observed at day -11. 

Our identification of an external layer of C in 
SN 1999aa remains tentative (based only on the contri- 
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bution at 6150 A), but seems plausible given the work 
of other authors on early spectra of SN 1991T-like and 
Branch-normal SNe. The presence of C would point 
to explosion hydrodynamics consistent with pure de- 
flagration of a Chandrasekhar mass WD (e.g. Nomoto, 
Thielemann, & Yokoi (1984)). 

Lentz et al. (2002) discussed the possibility that the 
absorption redward of Si n /16355 could be attributed 
to hydrogen mixed into the external layer from the WD 
companion. They modeled the presence of solar com- 
position material mixed in the unburned C+O layer of 
W7 to study the effect of different mixing depth. In 
some configurations they found that Ha could be eas- 
ily mistaken for C n /16580. We have explored this pos- 
sibility assuming the applicability of the resonance- 
scattering approximation, as illustrated in Fig. 21. 

The continua of the three models proposed have 
been shifted (adding an arbitrary constant) to match 
the data. The model on the top of Fig. 21 includes 
only Si n and, as noted in section 4.1, the minimum 
velocity is set by the photosphere velocity which is 
constrained by other ions to be v p h t= llxlO 3 km s _1 
and cannot be changed without adversely affecting the 
overall match. The model in the middle reproduces 
better the line profile thanks to the contribution of the 
Ha line. The parameters used for H are: v m i n =18xl0 3 
km s" 1 , v max =30xl0 3 km s _1 , t=0.08, T exc =15xl0 3 K 
and v e =5xl0 3 km s _1 . The optical depth used does not 
form any visible H/3 line around 4500 A, which is in 
agreement with our observations. Note, however, that 
for resonance scattering the net emission phenomenon 
is neglected and this could change the relative strength 
of Ha and H/3 (see e.g. Baron et al. (2000); Thomas 
et al. (2003)). The reproduction of the line profile at 
6150 A is indeed as good for Ha as for C n (bottom of 
Fig. 21). The identification of a given ion in a super- 
nova atmosphere based on the evidence of a single line 
is not definitive, therefore we consider the contribution 
from either of the two ions equally plausible. 

The presence of hydrogen in the external layers of 
the exploding WD would indicate that a single degen- 
erate progenitor system scenario is responsible for at 
least some of the observed SNe la. 

6.4. Carbon m or Hydrogen? 

In section 4.1 we included C in /14648.8 at photo- 
spheric velocities to match the line at 4500 A. The 
presence of unburned material deep into the IME at- 
mospheric region is seen in 3D deflagration models, 



where convective flows are shown to continuously re- 
fuel the inner layers with material that has original 
composition (Gamezo et al. 2003). Finding unburned 
material deep into the atmosphere is then theoretically 
possible. 

The first attempt to identify this line as C m was per- 
formed by Hatano et al. (2002) using the early spec- 
tra of SN 1997br, but they could not match the central 
wavelength of the line. They alternatively mentioned 
the possibility that He n could be responsible for the 
absorption. However, because of the high ionization 
energy of helium this could only be the case if there 
was a large degree of non-thermal ionization from the 
decay products of radioactive elements. 

Thomas et al. (2003) tentatively suggested the pos- 
sibility that the same line in SN 2000cx could be 
matched by a HV clump of H. In our spectrum at 
day 11 prior to maximum this would not be possi- 
ble because the optical depth necessary to match the 
line at 4500 A with H/3 would produce Ha in the 
Si ii /16355 A region, which is ruled out by the data 
as shown in Fig. 22 (SYNOW parameters for H: 
v min =24.5xl0 3 km s~\ v max =30xl0 3 km s _1 , r=2.5, 
T exc =15xl0 3 K and v c =3xl0 3 km s" 1 ). However, note 
that considering the net emission effect would weaken 
the Ha absorption feature. Furthermore, as pointed out 
in Thomas et al. (2003), because of the lack of a strong 
emission in the 6150 A region, the identification of 
the line at 4500 A as H/3 would be contingent upon 
a clumpy geometry and thus requires a 3D simulation 
to be fully explored. 

Because of the good line profile fit, C in remains 
our best candidate for the identification of the line at 
4500 A in the early spectra of SN 1999aa. Such a pos- 
sibility would strongly favor current 3D deflagration 
models of C+O WD since they produce this ion at all 
velocities, (Gamezo et al. 2003; Baron et al. 2003; 
Khokhlov 2000). Further confidence in this identifica- 
tion will require detailed radiative transfer modeling. 

6.5. High velocity Ca n 

The identification of high velocity burned material 
in supernova spectra can place constraints on the hy- 
drodynamics of the explosion. The presence of a HV 
component of Ca n is strongly supported by the blue 
components of Ca n IR triplet (Fig. 8). In the synthetic 
spectrum, the strength of the bluest of the two HV ab- 
sorptions features is almost twice that of the red one. 
This is not the case in the observed feature and could 
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be a sign of asymmetry in the high velocity ejecta. To 
reproduce the correct flux level, a three dimensional 
study of the geometry of the atmosphere (such as that 
computed for SN 200 lei in Kasen et al. (2003)) would 
be necessary. An optically thick HV clump viewed 
from a line of sight such that it just barely covers the 
photosphere can produce a weak IR triplet feature with 
two minima of equal depth (Kasen et al. 2003) such as 
seen in SN 1999aa. 

Evidence for a HV component of Ca n also comes 
from the parameters in our synthetic spectra needed to 
model observed Ca n H&K . To reproduce completely 
the broad profile of this line at day -11 and the blue 
minimum at -1 days the HV component is needed. 
When the temperature of the atmosphere drops, and 
the photosphere recedes into the inner layers, low ion- 
ization Fe-peak group lines become more important 
in the Ca H&K region. In this case, reproducing the 
line profile accurately becomes difficult. However, we 
could not match the Ca H&K absorption shown at day 
-11 as well as at day -1 by any means other than a 
HV Ca ii component. The evolution of this feature 
shows that the red component becomes stronger, and 
that by three weeks after maximum it dominates the 
entire feature. This would be consistent with a de- 
creasing optical depth of the HV Ca H&K as the at- 
mosphere expands. This is consistent with our inter- 
pretation, namely that the HV Ca n component could 
be formed in a clump. 

Recent spectropolarimetry measurements of type la 
supernovae (Howell et al. 2001; Wang et al. 2003; 
Kasen et al. 2003) have shown that, when present, the 
HV component of the Ca n IR triplet is usually po- 
larized. The study of the polarization parameters can 
yield important information on the geometry of the ex- 
plosion and thus can help in the development of 3D 
models. 

7. Conclusions 

We have presented new high signal to noise opti- 
cal spectroscopic data of SN 1999aa with good tempo- 
ral coverage between -11 and +58 days with respect 
to B-band maximum light. The overall evolution of 
the spectral features suggests an object with charac- 
teristics common to both SN 1991T-like and Branch- 
normal supernovae. 

By means of S YNOW synthetic spectra we have at- 
tempted to identify the absorption features in the spec- 
tra at days -11, -1, +5, +14 and +40 with respect to 



B-band maximum light. Highlights of this modeling 
include the presence of 

• C m at the photospheric velocity, 

• Possible C ii or H at high velocity 

• Ca ii IR triplet at high velocity, 

• doubly ionized Si, Ni and Fe, 

• confined IMEs and Fe n, 

• probable iron peak core at 10000 km s _1 . 

A schematic view of the resulting composition (in the 
velocity regime) is presented in Fig. 23 . 

The line identification for the earliest spectrum 
shows the presence of doubly ionized Si, Fe and Ni 
suggesting that the temperature in the outermost layer 
of SN 1999aa is higher than in normal SNe. After 
day 5, doubly ionized elements no longer form visi- 
ble absorption features, highlighting the cooling of the 
supernova atmosphere. The presence of IME above 
14000 km s" 1 excludes the possibility that SN 1999aa 
was the result of a sub-Chandrasekhar mass explosion. 

The broad absorption feature around 6150 A in 
the spectrum at day -11 can be matched compara- 
bly well by a weak component of Si n /I6355 plus 
HV C ii /16580 or Ha. The contribution of either C n 
or Ha helps in reproducing the velocity time evolu- 
tion of the absorption feature at 6150 A that otherwise 
(i.e. considering only Si n /16355) would remain un- 
explained. Definitive identification is not possible due 
to the lack of other absorption features produced either 
by C ii or H at the input optical depths. However, both 
the possibilities impose constraints either on the nature 
of the progenitor system or the hydrodynamics of the 
explosion. A hydrogen component would imply that 
SN 1999aa is the result of a single degenerate WD ex- 
plosion; a complete 3D calculation would be necessary 
to model it. Alternatively, a C n component can be re- 
produced by a pure ID deflagration model, but would 
not be consistent with delayed-detonation models. 

C m /14648.8 at photosphere velocity was able to 
reproduce the absorption around 4500 A, and no good 
alternatives were found for that feature. If this identifi- 
cation is confirmed by detailed radiative transfer mod- 
els, the use of three dimensional simulations of the ex- 
plosion would be necessary to describe the presence of 
unburned material down to the inner layers of the WD 
atmosphere. 
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At day -11 and at day -1, the profile of the absorp- 
tion at 3800 A requires a HV component of Ca to be 
reproduced in addition to the photosphere component. 
The HV component at day -1 also reproduces the ob- 
served minima on the blue side of the usual Ca IR 
triplet. For an accurate description of this feature, 
multi-dimensional simulations would be required. 

The spectral evolution of SN 1999aa in the photo- 
spheric phase shows that Fe n is confined below 15000 
km s _1 . IMEs populate a narrow velocity window 
above 10000 km s _1 . Similar evidence is found in 
other well known supernovae (SN 199 IT, SN 1997br 
and SN 2000cx) by studying the time evolution of 
the expansion velocity as computed from fits to the 
minima of Si n /16355. The comparison with other 
SN 1991T-like objects suggests that the transition be- 
tween IME to iron-peak dominant composition can oc- 
cur at slightly different phases. Tuning of the deflagra- 
tion to detonation transition density might reproduce 
this sequence. 

The analysis of the spectrum of day +40 indicates 
that the iron-peak core limit should be set to 10000 
km s -1 , similar to that of SN 199 IT (Fisher et al. 
1999). 

The origin of the differences between normal super- 
novae and SN 1991T-like or SN 1999aa-like objects 
depends on several factors. Our analysis of the optical 
spectra of SN 1999aa reveals that, among the present 
explosion models, none is able to reproduce each one 
of our findings. Higher temperature could account for 
some of the peculiarities (Nugent et al. 1995), but the 
evidence for high velocity components and unburned 
material at all velocities probably requires the develop- 
ment of full NLTE 3D explosion models. SN 1999aa 
was spectroscopically less extreme than other genuine 
SN 199 IT- like SNe, suggesting that perhaps a sin- 
gle model could eventually explain both normals and 
SN 1991T-like SNe as a continuous sequence. 
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Fig. 1.— SN 1999aa in its host galaxy NGC 2595. R.A. = 8/i27m42".03, Decl. = +21°29'14".8 (equinox 2000.0). 
B-band image obtained at NOT on 1999 February 13 UT with SN 1999aa indicated. The field is 6'.5 across. North is 
up, east to the left. 
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Fig. 2. — SN 1999aa spectral time sequence. The spectra are labeled in days relative to B-band maximum light. The 
© symbol marks telluric atmospheric absorptions. 
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Fig. 3.— The day -11 spectrum of SN 1999aa together with those of SN 1 99 IT, and normal SN 1990N andSN 1994D 
from Filippenko et al. (1992); Leibundgut et al. (1991); Patat et al. (1996). Each spectrum is labeled with the phase 
(days since B maximum). Line identification is explained in the text. The © symbol marks atmospheric absorptions. 
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Fig. 4. — Synthetic spectrum compared with SN 1999aa spectrum for day -11. SYNOW parameters are presented in 
Table 2. The region around 6150 A and 4500 A are enlarged in Figs. 5 and 6 to highlight the possible contributions of 
C ii and C ra. The mismatch of the continuum level above 5000 A has to be considered a limitation of the underlying 
black body assumption but does not affect the line identification. Ions responsible for features in the synthetic spectrum 
are labeled. 
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Fig. 5. — Synthetic spectra compared with SN 1999aa spectrum for day -11 around 6150 A. First model from the top: 
Solid lines, Ten = and data; Dashed lines, Ten = and Si n with v max = 30 x 10 3 km s _1 . Second model from the 
top: Solid lines, Ten + and data; Dashed lines, Ten + and Si n with v max = 30 x 10 3 km s _1 . The continuum level 
has been shifted to match the data. Ions responsible for features in the synthetic spectrum are labeled. The position of 
the minimum changes by about 50 A when the the C n component is included. The small notch around 6200 A is an 
artifact of the spectrum extraction and does not have to be considered a supernova feature. 
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Fig. 6. — Synthetic spectra compared with SN 1999aa spectrum for day -11 around 4500 A. Solid lines from top to 
bottom: rem = and data, rem + and data. 
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Fig. 7.— The -1 day spectrum of SN 1999aa together with those of SN 1991T, SN 1981B and SN 1994D from 
Jeffery et al. (1992); Branch et al. (1983); Patat et al. (1996). Each spectrum is labeled with the phase (days since B 
maximum). The © symbol marks atmospheric absorptions. 
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Fig. 8. — Synthetic spectrum compared with SN 1999aa spectrum for day -1. SYNOW parameters used are pre- 
sented in Table 3. The region around 3800 A and 8000 A are enlarged in Figs. 9 to highlight the contribution of 
Ca n HV. Above 6200 A the black body assumption fails in reproducing the correct flux values but it does not affect 
the identification of the lines. Ions responsible for features in the synthetic spectrum are labeled. 
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Fig. 9. — Synthetic spectrum compared with SN 1999aa spectrum for day -1 around 3800 A (left panel) and 8000 A 
(right panel). First model from the top: Solid lines, Tcaii hv = and data; Second model from the top: Solid lines, 
Tcaii hv + and data; The continuum level on the right panel has been shifted to match the data. Ions responsible for 
features in the synthetic spectrum are labeled. 
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Rest Frame Wavelength [A] 

Fig. 10.— The +5 day spectrum of SN 1999aa together with those of SN 1991T, SN 1990N and SN 1994D from 
Filippenko et al. (1992); Leibundgut et al. (1991); Patat et al. (1996). Each spectrum is labeled with the phase (days 
since B maximum). The © symbol marks atmospheric absorptions. 
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Rest Frame Wavelength [A] 

Fig. 11. — Synthetic spectrum compared with SN 1999aa spectrum for day +5. S YNOW parameters used are presented 
in Table 4. Above 6200 A the black body assumption fails in reproducing the correct flux values but it does not affect 
the identification of the lines. Ions responsible for features in the synthetic spectrum are labeled. 
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Fig. 12.— The +14 day spectrum of SN 1999aa together with those of SN 1991T, SN 1990N and SN 1994D from 
Phillips et al. (1992); Leibundgut et al. (1991); Mazzali et al. (1993). Each spectrum is labeled with the phase (days 
since B maximum). The ffi symbol marks atmospheric absorptions. 
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4000 5000 6000 7000 8000 9000 10000 

Rest Frame Wavelength [A] 

Fig. 13. — Synthetic spectrum compared with SN 1999aa spectrum for day +14. SYNOW parameters used are 
presented in Table 5. Above 6200 A the black body assumption fails in reproducing the correct flux values but it does 
not affect the identification of the lines. Ions responsible for features in the synthetic spectrum are labeled. 
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Fig. 14.— The +25 day spectrum of SN 1999aa together with those of SN 1991T, SN 1981B and SN 1994D from 
Phillips et al. (1992); Branch et al. (1983); Patat et al. (1996). Each spectrum is labeled with the phase (days since B 
maximum). The ffi symbol marks atmospheric absorptions. 
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Rest Frame Wavelength [A] 

Fig. 15.— The +40 day spectrum of SN 1999aa together with those of SN 1991T, SN 1981B and SN 1994D from 
Gomez, Lopez, & Sanchez (1996); Branch et al. (1983); Patat et al. (1996). Each spectrum is labeled with the phase 
(days since B maximum). The © symbol marks atmospheric absorptions. 
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Rest Frame Wavelength [A] 

Fig. 16. — Synthetic spectrum compared with SN 1999aa spectrum for day +40. SYNOW parameters used are 
presented in Table 6. A discontinuity at 10000 km s _1 has been introduced in the Co and Fe optical depths, indicating 
the possible iron-peak core limit. Ions responsible for features in the synthetic spectrum are labeled. 
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Fig. 17. — Doppler shift of the Ca H&K. The values for SN 1999aa are compared with those for other SNe la taken 
from Wang et al. (1996); Garnavich et al. (2001); Kirshner et al. (1993); Patat et al. (1996); Jha et al. (1999) and 
references therein. 
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Fig. 18. — Doppler shift of the Si n /J6355. The values for SN 1999aa are compared with those for other SNe la taken 
from Li et al. (1999, 2001a); Garnavich et al. (2001); Salvo et al. (2001) and references therein. The low value of 
SN 1999aa at day -1 1 is probably due to C n or H contamination. SN 2000cx, SN 1999aa and SN 199 IT maintain a 
constant velocity of 12000 km s _1 , 10100 km s _1 and 9400 km s _1 (dashed lines) respectively from maximum light to 
after day 20. 
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Fig. 19. — Doppler shift of the Fe m ,14404 . The values for SN 1999aa are compared with those for other SNe la 
taken from Li et al. (1999, 2001a) and references therein. 
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Fig. 20. — Doppler shift of the Fe m ,15129 . The values for SN 1999aa are compared with those for other SNe la 
taken from Li et al. (1999, 2001a) and references therein. 
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Fig. 21. — Synthetic spectra compared with SN 1999aa spectrum for day -11 around 6150 A. First model from the 
top: Solid lines, t h = 0, Ten = and data. Second model from the top: Solid lines, t h + 0, Ten = and data. Third 
model from the top: Solid lines, Ten + 0,t h = and data. The continuum level has been shifted to match the data. 
SYNOW parameters for C n are presented in table 2 and those for H are presented in the text. Ions responsible for 
features in the synthetic spectrum are labeled. 
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Fig. 22. — Synthetic spectra compared with SN 1999aa spectrum for day -11 around 4500 A (left panel) and 6150 A 
(right panel). First model from the top: Solid lines, t h + 0, rem = and data. Second model from the top: Solid lines, 
th = 0, rem + (as second model in Fig. 6) and data. In the resonance scattering approximation the optical depth 
necessary to reproduce the observed absorption at 4500 A with H/3 produces a too strong Ha that could be weaken 
considering the net emission effect. S YNOW parameters for C m are presented in table 2 and those for H are presented 
in the text. Ions responsible for features in the synthetic spectrum are labeled. 
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Fig. 23. — Schematic representation of the composition of SN 1999aa in velocity space as inferred from the SYNOW 
synthetic spectra of day -11, -1, +5, +14 and +40. Velocities increase radially in units of 10 3 km s _1 . The velocity 
ranges (v m! „-v max ) for 12 ion species used in the synthetic spectra are marked with solid lines inside the circle section 
labeled with the ion name. High velocity components are marked with dashed lines. In each ion section the five epochs 
analyzed are shown. For each epoch the extent of the photosphere is marketed with a thin solid line. The contributions 
of C n, C ra, [O n] and H are not shown. For details see section 4. 
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Table 1: Data set specifications. 



JD 


Epoch' 


Telescope 


Instrument 


A Range 1 * 




(S/N) b 


Comments 


-2400000 


ref B max 






[A] 


[A] 






51223.38 


-11 


NOT 2.6m 


ALFOSC 


3503-9295 


6.1 


52 


5998 c 


51227.84 


-7 


APO 3.5m 


DIS 


3686-10268 


6.8 


75 


5490 c 


51231.35 


-3 


NOT 2.6m 


ALFOSC 


3508-9692 


6.1 


41 


6097 c 


51234.71 


-1 


Lick 3 m 


KAST 


3478-9975 


3.1 


43 


5430 c 


51239.68 


+5 


MDM 2.4m 


MARK III 


3872-8807 


5.3 


45 


e 


51239.68 


+6 


MDM 2.4m 


MARK III 


3871-8808 


5.3 


61 


e 


51247.63 


+ 14 


APO 3.5m 


DIS 


3628-10333 


6.8 


86 


5887 c 


51253.37 


+ 19 


NOT 2.6m 


ALFOSC 


3375-9994 


6.0 


63 


6057 c 


51258.51 


+25 


CTIO 4m 


RCSP 


3320-9208 


2.0 


40 


d 


51261.54 


+28 


CTIO 4m 


RCSP 


3374-9232 


2.1 


42 


d 


51266.53 


+33 


CTIO 4m 


RCSP 


3263-9199 


2.1 


42 


d 


51273.62 


+40 


APO 3.5m 


DIS 


3641-10115 


6.7 


42 


5877 c 


51282.67 


+47 


Lick 3 m 


KAST 


3512-9940 


3.1 


18 


5462 c 


51286.63 


+51 


APO 3.5m 


DIS 


3621-10166 


6.7 


68 


5874 c 


51293.71 


+58 


Lick 3 m 


KAST 


3705-7963 


2.1 


11 


5390 c 



"Average wavelength-bin size. 

6 Average signal-to-noise ratio per wavelength bin. 

beginning of red channel, [A]. 

^Negligible 2 nd order contamination. 

e Possible 2 nd order contamination above 7500 A. 

f In Rest Frame. 

fWhere S/N > 5. 
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Table 2: SYNOW parameters for day -11. The fit is shown in Fig. 4. v p h ot =l 1000 km s -1 , r bb = 13700 K. (HV) next 
to Ca n stands for high velocity component. 



Ion 


T 


Vmin 

lO^ms- 1 


Vmax 

lO^ms- 1 


Texc 
10 3 K 


V e 

10 3 km s- 1 


Ca ii 


1.15 




19.5 


15 


5 


Ca ii (hv) 


1.45 


19.5 


30 


15 


10 


Sin 


0.1 




15 


15 


5 


Sim 


0.3 




16 


15 


5 


Con 


0.25 


14 


30 


10 


5 


Fern 


0.63 




30 


10 


5 


Ni in 


2 


14 


30 


10 


5 


Cn 


0.01 


19 


30 


15 


5 


Cm 


0.5 




13.5 


15 


5 
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Table 3: SYNOW parameters for day -1. The fit is shown in Fig. 8. v phot =9500 km s -1 , r bb = 13500 K. (HV) next to 
Ca ii stands for high velocity component. 



Ion 


T 




Vmax 


7"exc 


V e 






10 3 kms-' 


10 3 kms-' 


10 3 K 


10 3 km s _1 


Ca ii 


3 




18 


8 


5 


Ca ii (hv) 


2.5 


22 


25 


10 


5 


Oi 


0.2 




30 


8 


5 


Sin 


0.9 




15 


8 


5 


Sim 


0.25 




15 


8 


5 


Sn 


0.65 


10 


20 


8 


5 


Mg ii 


0.1 


10 


20 


8 


5 


Fe ii 


1 


10 


15 


8 


5 


Co ii 


0.3 


13 


30 


8 


5 


Nai 


0.25 




30 


8 


5 


Fe in 


0.8 




30 


8 


3 


Nin 


0.1 




30 


8 


5 
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Table 4: SYNOW parameters for day +5. The fit is shown in Fig. 11. v phot =9000 km s -1 , r bb = 12000 K. 



Ion 


T 


Vmin 

lO^ms- 1 


Vmax 

10 3 km s- 1 


10 3 K 


V e 

10 3 km s- 1 


Can 


20 




30 


5 


3 


Oi 


0.1 




30 


5 


3 


Si ii 


1.5 


10 


15 


5 


3 


Sim 


0.05 


10 


15 


5 


3 


Sn 


0.45 


10 


20 


5 


3 


Mg ii 


0.6 


10 


20 


5 


3 


Fe ii 


1.3 


11 


14 


5 


3 


Co ii 


0.8 


11.5 


30 


5 


3 


Nai 


0.4 




30 


5 


5 


Fe in 


0.09 




30 


5 


3 


Ni ii 


0.2 




30 


5 


3 
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Table 5: SYNOW parameters for day +14. The fit is shown in Fig. 13. v phot =8500 km s" 1 , T bb = 10500 K. 



Ion 


T 


lO^ms- 1 


Vmax 

10 3 km s-> 


T 

1 exc 

10 3 K 


v e 

10 3 km s" 1 


Can 


20 




30 


5 


3 


Oi 


0.1 




30 


5 


3 


Sin 


1.5 


10 


14 


5 


3 


Fe ii 


4 


9.5 


13.5 


5 


3 


Con 


1.7 


10 


30 


5 


3 


Nai 


1.2 




30 


5 


5.5 


Nin 


5 




10 


5 


3 
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Table 6: SYNOW parameters for day +40. The fit is shown in Fig. 4. v phot =8000 km s -1 , r bb = 8000 K. 



Ion 


T 


Vmin 

lO^ms- 1 


^max 

10 3 km s- 1 


T 

* exc 

10 3 K 


v e 

10 3 km s- 1 


Ca ii 


300 




18 


8 


3 


On 


1 




30 


10 


10 


Fe ii 


160 




10 


5 


3 


Fe ii 


0.2 


10 


13.5 


5 


3 


Co ii 


60 




10 


5 


3 


Co ii 


0.2 


10 


30 


5 


3 


Nai 


0.9 




30 


5 


5 
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